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ABSTRACT 

Context. A large sample of T Tauri stars exhibits optical jets, approximately half of which rotate slowly, only at ten per 
cent of their breakup velocity. The disk- locking mechanism has been shown to be inefficient to explain this observational 
fact. 

Aims. We show that low mass accreting T Tauri stars may have a strong stellar jet component that can effectively brake 
the star to the observed rotation speed. 

Methods. By means of a nonlinear separation of the variables in the full set of the MHD equations we construct semi- 
analytical solutions describing the dynamics and topology of the stellar component of the jet that emerges from the 
corona of the star. 

Results. We analyze two typical solutions with the same mass loss rate but different magnetic lever arms and jet radii. 
The first solution with a long lever arm and a wide jet radius effectively brakes the star and can be applied to the 
visible jets of T Tauri stars, such as RY Tau. The second solution with a shorter lever arm and a very narrow jet radius 
may explain why similar stars, either Weak line T Tauri Stars (WTTS) or Classical T Tauri Stars (CTTS) do not all 
have visible jets. For instance, RY Tau itself seems to have different phases that probably depend on the activity of the 
star. 

Conclusions. First, stellar jets seem to be able to brake pre-main sequence stars with a low mass accreting rate. Second, 
jets may be visible only part time owing to changes in their boundary conditions. We also suggest a possible scenario 
for explaining the dichotomy between CTTS and WTTS, which rotate faster and do not have visible jets. 

Key words. MHD, Stars: pre-main sequence. Stars: winds, outflows. Stars: mass-loss. Stars: rotation, ISM: jets and 
outflows 

1. Introduction breakup speed (|Matt et all . |2010[ ). This indicates that a 

' . very efficient mechanism is at work to remove the angular 

7^. ; 1.1. The angular momentum problem, observational facts momentum in these stars; the nature of this mechanism is 

. !^ ' Observations of star-forming regions show that several T controversial. 

^ . Tauri stars are associated with well coUimated jets. These iKundurthv et all (|2006D . see also lEdwards et al.l (|1993[ ). 
H ; jets are made of stro ngly accelera ted plasma that travels at claimed that CTTS are slow rotators with circumstellar 
. P. , a few hundred km/s JBaliyl, 120091) and is niagneti cally colli- disks, while WTTS are fast rotators without accretion 
mated at large distances (Dougados^LaD, |2000|). The jets disks. However, the association of CTTS with slow rota- 
are usually associated with Classical T Tauri Stars (CTTS) , tors and WTTS with fast ones is not a s yet observa tion- 
which are low mass stars < 2Mq in their late stages of pre- ally confirmed. Mo reover, IStassun et all (|1999l |2001[ ) and 
main sequence evolution, showing strong evidence of the iRebull et all (|2004f l could not find a strong correlation of 
presence of a surrounding accretion disk. The outflow is fast rotators with stars that have d ispersed their disks. 
usually thought t o be e jected from the Keplerian disk (e.g. Conv ersely, new results (iRebull et all . [2006; Hcrbst et all 
ICabrit fc A^drj |l99lD. However, a significant part of the |2007l : ICieza fc Ballbe^ . I2007D founda bimodal distribution 



jet, if not all, may be ejected by the star itself, at least for where stars with disks rotate slower than stars without 
the lower mass accreting T Tauri stars. disks. However, the two distributions clearly overlap. The 
Observations also reveal that approximately half of the prevailing theory has been that magnetic interaction be- 
T Tauri stars rotate slowly, at about 10% or less of their tween the star and its circumstellar disk regulates the stel- 
lar rotation periods. Attempts to observationally confirm 



Send offprint requests to: christophe.sauty@obspm.fr this conjecture have produced mixed results. Hence, from 
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the observational point of view, the presence of an accretion 
disk may not be enough to explain stellar angular momen- 
tum removal. 

In order to understand slow stellar rotators. ISchatzmanI 
(|l962ll suggested that the magnetic braking of a stellar wind 
would be sufficient even with a low mass l oss rate, and this 
has been explored by several authors (e.g. IWeber fc DavisI , 
[T967. Mcstel , 1968a, 1968b) . Following this track, it is es- 
sential to address stellar wind models in detail, with and 
without the presence of an accretion disk, to see how ef- 
ficient the wind is in extracting angular momentum from 
the star. Moreover, recent observations confir m that stellar 
winds are present in at least 60% of CTTS (|Kwan et all . 
20071). as was alread y suggested by various authors (e.g. 



Edwards et aD,'2003'). 



[Decampli (1981) showed that t/ierTna/ stellar winds can- 
not support jets with a high mass loss rate, simply because 
this would require too high temperatures to be physical. For 
low mass loss rates, less than lO^^Mo/yr, he showed that a 
thermally driven or line-driven radiative wind can support 
the formation of the jet. For a higher mass loss rate, such 
a driving is still possible provided that there is an extra 
supply of pressure from Alfven waves. If the temperature 
of the wind is around one million degrees, the pressure can 
drive the jet up to 300 km/s. However, the pressure may 
exceed the thermal pressure by a large amount if it includes 
ram pressure or turbulent magnetic pressure. The situation 
is similar to the well known case of the solar wind, where to 
obtain wind speeds of several hundred km/s with a temper- 
ature of a few million degrees one requires extra pressure 
from either turbulence or kinetic effects. These effects yield 
an effective t empe rature of around ten million degrees (see 
lAibeo et al.l . l2007l and references therein). 

Ignoring the important question regarding t he nature 
of the pressure and following the early study of iDecamplil 
(|198lD . most authors disregarded stellar wind models up 
to now. This explains why disk wind models have been 
extensively studied. 

1.2. Disk wind theory 

In disk wind models, the main source of acceleration is 
the m agnet o-centrif ugal driving, as proposed originally by 
iBlandford fc Pavnd (l982 ). who presented the first MHD 
radial self-similar so lutions extending t h e earl ier hydrody- 
namical solutions of iBardeen fc Bergeil (|1978[ ) for galactic 
winds. Since then, the ir original model has been improve d 
in various ways (e.g., Li et al., 1992; Contopoulosl |1994| 



IVlahakis fc Tsineanos. 1998: Vlahakis et al. . 20001 ) . In par- 
ticular iFerreiral (|l997t ) has shown that a consistent treat- 
ment of the connection with the disk is possible. From the 
theoretical point of view however, self-similar disk wind 
models are not consistent ly describing the in ner part of 
the jet close to its axis (jGracia et al.l . |2006[) . The solu- 
tion has to be cut at the inner and outer edges of the 
disk. From the observational point of view, c old disk winds 
turned out to be too rapid and too light (iGarcia et al.l . 
l2001btl . This last point can be overcome with warm disk 
solutions and explains the high speed, dense powerful jets 
of early TTauri stars like DG Tau, as well as their rotation 
([Ferreira et al.l . l2006t ). DG Tau is probably one of the most 
powerful sources with an accretion rate of 10~^Afo/yr and 
a corr esponding mass loss r ate that is about 10% of this 
value (jHartigan et al.l . ll995l) . In this case, disk wind solu- 



tions seem more adapted. Recent observations suggest that 
DG Tau is variable on short timescalcs with an accretion 
rate dropping down to 10~^Mo/yr (Beck ct al., 2010). 

Simultaneously, several numerical simulations have 
investi gated the jet l a unchin g of these disk winds 



e.g.. lOuved fc Pudritj. Il997t lUstvugova et all . 120001: 



iKras nopolsk v et al.l . l200l ICasse fc Keppensl . 120041 ). By in 



vestigating various boundary conditions and initial configu- 
rations, they confirmed that disk winds could indeed be ac- 
celerated and collimated by their own m agnetic field. More 
recently several authors (iMatsakos et al .. 2008; S tute et all 
120081 : IMatsakos et al.l . '2009) have shown that analytical so- 
lutions such as the radially self-similar ones are structurally 
stable. This result indicates that these solutions are good 
complementary tools to numerical simulations. 

According to theoretical arguments 



(iPudritz fc Norman . 1986) and numerical simulations 
(e.g. iMeliani et al.l . 12006. ) . disk winds can remove most of 



the angular momentum of the accreting plasma. However, 
disk winds need the presence of an inner stellar jet 
component to be globally consistent. 

From the theoretical point of view, it was initially 
proposed that fieldlines of the stellar magnetosphere an- 
chored into the disk could efficiently brake the central 
star, in the so-called disk-locking mechanism proposed by 
Choi fc H erbst (1996). In principle a fieldline rooted into 
the star will slow down (respectively speed up) the central 
star rotation, if it connects to a region of the disk rotating 
at lower rotational speed (respectively higher speed) after 
(respectively before) the corotation radius. Using a simple 
analyti c al mod el based on the Ghosh fc Lamb mechanism 
(|1979al. Il979b[ ). iMatt fc PudTit^ (|2005f ) showed however 
that the disk-locking efficiency is reduced because it leads to 
an opening of the fieldlines that disconnects the disk from 
the star. Instead they proposed that angular momentum 
can be lost via accretion-powered stellar winds provided 
the mass loss rate is 1 to 10 % of the accretio n rate as usu- 
ally inferred from observations ([Calvetl . 11998 ) . Their results 
were confirmed via numer i cal simulations (lMa.tt fc PudritzL 
I2008al 120085 IMatt et al] . l2010l ). lKiiker et al.l (l2003l) . using 
similar simulations, concluded that even when the disk- 
magnetosphere connection is not disrupted, the magnetic 
field could slow down the central object but the domi- 
nant disk torque still spins up the star. As mentionned by 
JRomanova et al.l (|2009[ ). these conclusions are not defini- 
tive because long-term simulations are required to establish 
the validity of this result. Ultimat ely the authors insist on 
the essential role of the axial jet. 'Zanni fc Ferreira! (|2009l) 
performed various additional simulations where the disk- 
locking is not disrupted by the building of large toroidal 
field. However, even in this situation, the magnetic accre- 
tion torque is not sufficient to cause an efficient spin-down 
of the star. 



1.3. Stellar Jet theory 

As a conclusion of the precedent subsection, the presence 
of a pure stellar wind in addition to a possible outer disk 
wind is required to effectively remove the excess of stellar 
angular momentum. 

Besides, more evolved T Tauri stars like RY Tau seem 
to have weaker jets with accretion rates on the order 
of only a few 10~^MQ/yr, which may be an indication 
that disk winds are not always absolutely necessary for 
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the weaker and more evolved jets. We see that with a 
mass loss rate of only 10% of this value, these jets are 
more difhc ult to analyze. From spectroscopic observa tions 
of RY Tau, iGomez de Castro fc Verdugol (|20niL I2007D sus- 
pected the presence of a small-scale pure stellar jet be- 
cause the UV emission lines originate in a region that 
is too sn iall to be produced by a disk win d. New obser- 
vatio ns (|St-Onge fc Bastiem 120081 : lAgra-Amboage et all . 
I2009D clearly show a microjet in this faint object. It is in- 
teresting to investigate if these (micro- )jets can be modeled 
through stellar winds or a combination of a stellar wind and 
a sub Keplerian disk wind. 

This combination of a two -component outflo w is un- 
der n umerical investiga ti on (iMatt et al.l 120031 : iKoidd . 



120031: iKiiker et al.l. 120031: 



simulations (IMeliani et al. 



Ma tsakos et al.'. '2008V Recent 
2006; Romanova et al., 2009; 



iMatsakos et al.l . l2009t iFendti . l2009i ) also show the interplay 
of mixing between the two components. Sever al a uthors, 
e.g.lKukeret al.l (|2003[ ) . iRomanova et al.l (|2009D and lFendtl 
(|2009t ) , have shown in their numerical simulations that the 
disk-magnetosphere connection produces strong intermit- 
tent outflows, the results of flares and reconnection. This 
is somewhat similar to coronal mass ejections in the so- 
lar wind and in microquasars. Though important to model, 
those outflows are not necessarily more important on the 
long-term duration compared to the continuous underly - 
ing stead y flow, as suggested by iRomanova et al.l (J2009t) . 
Moreover. IMatsakos et al] ()2009t ) have shown that the vari- 
ability of the source does not destroy the continuous steady 
jet. Indeed, they use the first solution presented in this pa- 
per, as initial condition for the inner spine jet of their simu- 
lations, but use a polytropic equatio n of state. In their pre- 
vious paper (jMatsakos et al.l . l2008l ) they have shown that 
changing from non polytropic to polytropic would reduce 
the size of the radius of the jet, without drastically affecting 
the outflow behavior, however. 

However, numerical simulations are usually time- 
consuming to such an extent that they cannot explore a 
wide range of parameters. Moreover, they do not really 
start ejecting mass directly from the star but rather at a 
given height between the sonic and the Alfven surfaces. 
Therefore, it is of absolute interest to model the stellar jet 
that originates in the central star. 

All these problems can be overcome by studying 
MHD solutions obtained via a nonlinear separation of 
the variables. These so-called self-similar solutions con- 
tain as special cases all known MHD outflow solutions 
(e.g., IVlahakis fc: Tsinganosl Sautv et al.. 2002, hereafter 
STT02). As shown in lSautv fc Tsin~os (19941 . hereafter 
ST94, these self-similar models can be seen as a combi- 
nation of stellar wind and X-wind models. Convers ely to 
formal X-winds (IShu et al.l . [19941: IShang et al.l . l2002D . they 
do not require that all fieldlines emerge from a single point 
of the disk in the form of a fan. Instead there is a smooth 
transition between the stellar wind and t he disk jet. They 
are also shown to be structurally stable ([Matsakos et al.l . 
[2008; 2009). This was far from obvious because the heat- 
ing function is not polytropic as usua lly in most a nalyt- 



ical solutions. In a series of papers - Sautv et al.l (|1999[ ) 
hereafter STT99, STT02, and lSautv et al.l (|2004D . hereafter 
STT04 - we have systematically explored the full parame- 
ter space of this problem. In ST94 we have shown that the 
double disk and star component was promising. However, 
the solutions presented there had mass loss rates that were 



too low. Here we present solutions where the stellar jet 
is both under-dense and under-pressured compared to the 
disk wind in the launching region. These solutions have the 
advantage to adapt themselves not only to CTTS, which are 
connected with the surrounding accretion disk, but also to 
WTTS where there is either no accretion disk or no con- 
necti on with the disk itself (Walt er et al.[ 119881 : iBertoutl 
I1989D . Although winds from WTTS are diflficult to mea- 
sure, mass loss rat es between 10~^^ an d 10~^''MQ/yr have 
been obtained by lAndre et al.l (|1992D . Winds of WTTS 
mig ht be similar to the solar wind in terms of collimation 
fsee lAibeo et al.[ l2007l and references therein). However, 
these young stars have higher mass loss rates and rota- 
tional speeds than our sun and might as well have strongly 
collimated jets independently of the presence of an accre- 
tion disk. Their rotation periods range from 0.6 to 24 days 
with peaks near two and eight days (jHerbst et al.[ |2007() 
while the solar rotational period is 24.5 day. Thus, WTTS 
may well produce self-collimated jets that nevertheless are 
so weak that they cannot be detected. 

Wc will use the ST94 meridionally self-similar model to 
examine the contribution of the stellar component to the 
overall jet and how efHciently it brakes the star. We first 
explore how the model parameters can be constrained from 
the observations of CTTS jets (Sect. |2|). Then, we focus 
our attention on two different classes of solutions (Sections 
m and [SJ and explore the efficiency of these winds to spin 
down the central object. In the final section (Sect. we 
discuss how these two classes of solutions may explain the 
dichotomy between CTTS and WTTS and how they can 
be adapted to also model multiphases of T Tauri faint jets, 
such as in RY Tau. 



2. Meridional self-similar solutions 

We summarize the main assumptions of our meridion- 
ally (0—) self-similar treatment of the MHD equations in 
Appendix [SI . More details can be found in STT94, STT99 
and STT02. 



2.1. Parameters and variables 

Using spherical coordinates {r,9,ip), all quantities are nor- 
malized to the Alfven surface along the rotation axis, which 
is by assumption a sphere of radius r = r^. The dimension- 
less radial distance is denoted hy R — rjr^, while B,, K 
and p, are the poloidal magnetic field, velocity and den- 
sity along the polar axis at the Alfven radius r*, with 
K' = S2/47rp,. 

The original system of MHD equations reduces to two 
coupled partial differential equations for the density and the 
magnetic flux. Then, the momentum conservation law pro- 
vides three ordinary differential equations, which together 
with Eq. (|A.12p can be solved for the four variables ]VP{R), 
F{R), H(i?) and G{R), which are the square of the poloidal 
Alfven Mach number, the expansion factor, the dimension- 
less pressure, and the dimensionless cross section radius, 
respectively. 

The model is controlled by the following four parame- 
ters. 

— The parameter d governs the non spherically symmetric 
distribution of the density. For 6 > {S < 0) the density 
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increases (decreases) by moving away from the polar 
axis. 

The parameter k controls the non spherically symmetric 
distribution of the pressure. For k > (k < 0) the gas 
pressure increases (decreases) by moving away from the 
polar axis. 

The parameter A is related to the rotation of the 
poloidal streamlines at i? = 1 (i.e., at the Alfven sur- 
face) . 

The parameter v is defined by the gravitational field 
whose acceleration can be written as 



GM 



IV^ 



i?2 



(1) 



where A^ is the central gravitating mass. This extra 
parameter v is the ratio of the escape velocity to the 
poloidal flow velocity on the polar axis at _R = 1 (i.e., 
at the Alfven surface), 



2gM 



2.2. Collimation efficiency of the magnetic rotator 



(2) 



By integrating the momentum equations along a fieldline, 
we obtain the conserved total energy fiux density per unit 
of mass flux density. This is equal to the sum of the kinetic 
and gravitational energies together with the enthalpy and 
net heating along a specific streamline. In the framework of 
the present meridionally self-similar model, the variation 
of the energy from one line to the other gives an extra 
parameter (STT99): 



Ar 



[GRf 



El 

4 



1 






A2 fAP-G^ 



G2 V 1 - M^ 



2X' 



R 

1 - Af 2 



(3) 



which is a constant on all streamlines (ST94). 

Physically, e is related to the variation across the field- 
lines of the specific energy that is left available to colli- 
mate the outflow once the thermal content converted into 
kinetic energy and into balancing gravity has been sub- 
tracted (STT99). 

We can express e/2\^ in terms of the conditions at the 
source boundary Tq (see STT99 for details). 



e _ ^Poynt.,o + £^R,o + A£^q 

2A2 Emr 



(4) 



where Emk is the energy of the magnetic rotator (see Eq. 
2.5a in STT99), i5poynt.,o is the Poynting energy and E-r^o 
is the rotational energy at the base. Aii^Q is the excess or 
deficit of the gravitational energy (per unit mass) that is 
not compensated by the thermal driving, on a nonpolar 
streamline compared to the polar one. 



AS? 



GM 



1 



To{a) 



To (pole) 



GM{5 



Kja 



1 + Sa 



(5) 



For e > collimation is mainly provided by magnetic 
terms, while for e < the outflow can be conflned only by 
thermal pressure, which is not possible for overpressured 



flows. Accordingly, in STT99 we classified fiows with posi- 
tive or negative e as efRcient or inefficient magnetic ro- 
tators, respectively (EMR or IMR). 

Considering the fairly large opening of stellar jets, we 
expect to find values of e close to zero. We have shown in 
STT99 that this corresponds to large asymptotic jet radius 
of magnetically confined outflows. 

3. Constraining the model parameters from the 
observations 

3.1. Stellar constraints 

Though their jets are fainter than those of class I and early 
CTTS, low mass accreting classical T Tauri stars exhibit 
jets as many other class II objects. Th ey have well measured 



mass, radius and rotation frequ ency (Hartigan c t al.l. 11995 : 



Bouyier et a!] . 1X9971: iHerbst fc Mundt . .2005: .Marilli et al 
2003). We suspected that one good such example of a 



faint jet in an evol ved CTTS would be, am ong oth- 
ers, the RY Tau jet (|G6mez de Castro fc Verdu eo. 2001j) 
This has been recently observed and extensively stud 



ied (|St-Onge fc Bastienl[2008tlAgra-Amboage et alll2009[ ). 
Classical T Tauri stars are low mass stars between 0.5 and 
2 solar masses, with typical radius ranging from 2 to 3 so- 
lar radii (Bertout, 1989; Mora et al., 2001). They are typ- 
ical precursors of solar-type stars. Rotational speeds are 
usually assumed to be around one tenth of their breakup 
speed, namely typically around 15 kms~^ for CTT S and 
around 20 kms"! for WTTS (|Bouvier et all . [1993), but 
with large variations and uncertainties. In particular, there 
is a large discrepancy between spectroscopic measurements 
of projected rotational velocities and periods measured us- 
ing photometry. 

3.2. Jet constraints 

For some CTTS we have measurements of various physi 
cal qua ntities in the optical jet itself tHartigan e t.al 
199" 



Lavallev-Fouauet et al. . 19991: Cabrit 



iSt-Onge fc Bastienl . 120081 : lAgra-Amboage et al 



120071: 
I2009h . 

Observations give mass loss rates from 10~^ down to 
lO~'^'^MQ/yr. On the other hand, low mass accreting 
stars are also associated with low mass loss rates, close 
to the minimum value. The averaged asymptotic speeds 
are on the order of 100 km s~^ to 300 kms~^ (jBallvl 
,2009). Microjets are hardly resolved transversely. Thus, 
this speed might be a mixing of the high- and low- velocity 
components sim i larly to what is observed in DC Tau 
(jAnderson et al.l . l2003f l. For the high- velocity component 
we can easily take a higher value and assume 400 km s~^ 
along the polar axis, where we expect to find this maximum 
value. This remains of course arbitrary but is consistent 
with the model developed here. In that respect, more 
precise measurements from observat ions are also needed t o 
assess the velocity profile in the jet (jGiinther et al.l . l2009t) . 
We assume that the jet carries a typical mass 
loss rate of a few tim es 10~^ M ^/yr. Following 
iGomez de Castro fc Verdtigol (|200lD . we took 

3.1 X lO~^M0/yr for the mass loss rate of RY Tau. 
Electronic densities can be directly measured only in 
shocks. Typically they are around 10* cm~'^. However, 
the inferred electronic density is not unambiguously 
determined. 
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First the jet radius may be much smaller than the 
shock radiu s, as shown in many numerical simulations (e.g. 
Matsakos et al. . 20091) . As mentioned in the introduction, 
Matsakos et al.l ( 20091 ) used as initial condition our solution 
for the inner spine jet but with different thermodynamics. 
Switching to a polytropic equation of state reduces the ra- 
dius of this inne r part of the ou tflow by a significant fac- 
tor, as shown in iMatsakos et all ( [2008 ) . Even with such a 
small radius of the spine jet, the long term variability of the 
stellar mass loss rate produces knots ten times larger than 
the jet radius, typically around 5 to 10 AU. We expect that 
keeping the turbulent heating along the flow would produce 
even larger knots. 

Second, not all electrons are accelerated to sufficiently 
high energies. Large uncertainties may also arise because 
we do not precisely know the ionization fraction of the jet. 
Thus, densities and temperatures cannot be calculated pre- 
cisely. The solutions presented in this paper can give a fairly 
accurate description of the outflow dynamics. In order to 
make a full synthetic emission map, further work is nec- 
essary starting from the dynamics given from these solu- 
tions, calculating the amount of mechanical heating and 
self-consistently deducing the emission lines by tuning the 
ionization fra ction and the kinetic temper ature. This has 
been done by iGarcia et al.l ()2001al . l2001bf ) for disk winds 
and will be similarly done for our solutions in a forthcom- 
ing study. 

Finally, to have efficient magnetic braking, we looked 
for MHD solutions with long magnetic lever arms. We as- 
sume that a typical lever arm, the Alfven cylindrical radius, 
is around 10 times the stellar radius, uia = IQvjo- For the 
meridionally self-similar model that we use this gives ap- 
proximately r* = 10 To or Rq = 0.1. 

As shown in lSpruid ()1997[ ). if the accreting mass rate is 
Mace, the angular momentum is accreted at a rate Jacc = 
-Macc'H^o^- Similarly, the angular momentum loss rate in the 
wind is Jwind = Afwind^^a^, bccause the plasma corotates 
up to the Alfven cylindrical radius Wa ■ If the entire angular 
momentum from the disk is removed by the wind, we can 
equal these quantities Jwind = Ja.cc- Combining the above 
we get 



Ma, 



Wa 



(6) 



If uja = 10 vjo, we deduce from the above equation that the 
mass loss rate is only 1% of the mass accreting rate, which is 
a fairly low value. The wind is probably more efficient than 
this 1% at extracting mass from the disk. This means that 
even a relatively small wind outflow can remove angular 
momentum from the disk very efficiently. 

Alternatively, the wind can also remove angular momen- 
tum from the star if the magnetic fieldlines and streamlines 
are anchored in the stellar photosphere instead of the disk. 
In that case, we can calculate the braking time to remove 
the stellar angular momentum as 



Js\ 



T ■■ 



•^wind 

where 

Jstar = kflMstm-ro, 



(7) 



(8) 



and 



Jvj 



2 V'out 



(9) 



with k the dimensionless inertial const ant of the star (?). 
For E g. © we used Eq. (3.19b) of ISautv fc TsinganosI 
(IT99I . 

We use k = 0.06 as for the Sun 
([Schwartz fc Schubert et~all . I1969D wher e 90% of its 
mass is within 50% of its radius ()Coa 20001) . For a T Ta uri 
star with a bigger convective zone (jBouvier et al.l . 119971 ) . k 
may be larger, up to 0.2. This maximum value is obtained 
for a fully convective star . ^_^ 

Substituting ft from Sautv fc Tsinganoi (|1994D in Eq. 
^ , we obtain the following expression for the braking time 
of the star 



2/cMstarrg 



V'o 



IM, 



wind 



Ckout^Vl + Sao 



with 



-00 



Sk' 



Sa„ 



\3/2 



!]• 



(10) 



(11) 



For instance, we know that CTTS have lifetimes on the 
order of a million years, while WTTS last on this phase 
around ten million years. The magnetic braking should ac- 
cordingly be efficient in removing most of the accreted and 
stellar angular momentum during this period, independent 
of the uncertainty factor we have on. 

3.3. Stellar parameters 

To summarize, we will use a mass loss rate, stellar mass, 
radius and rotation given respectively by 

- Mwind = 3.1x lO-^Mo/yr, 

- M ^ 1.5 Mq , 

- ro = 2.4rQ, 

- n = 5.99 10^6 rad/s or V^^o = 10 km s^^ . 

These values are typical of CTTS and close to the ones of 
RY Tau that we will study in more detail, except for the 
rotation rate, which is much lower than the usual value. 
We argued, however, that the rotation rate of RY Tau is 
subject to caution because of the obscuration of the star by 
the disk, as discussed below. Instead, we prefer to keep a 
value close to that of most CTTS and remain more general. 

Taking these constraints into consideration, we explored 
the parameter space and constructed various solutions us- 
ing the procedure explained in STT02. Among the solutions 
adapted to T Tauri jets, we selected two different sets of 
model parameters that fit the RY Tau case particularly 
well. The first solution has a longer lever arm, correspond- 
ing to our initial guess, and a large opening angle for the 
jet. This should be close to the observed values for CTTS 
with a j et, such as the jet o f RY T au observed in February 
2005 bv lSt-Oneefc BastienI (|2008l) . In this solution, the jet 
is magnetically coUimated into a cylindrical shape without 
any oscillation close to the star. 

By a slight change of the pressure, we were able to get 
a second solution that is completely different from the pre- 
vious one. Of course, a change in one of the parameters 
implies a search for a new set of the other parameters to 
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keep the same observational constraints. In this second so- 
lution, the lever arm is reduced, as is the opening angle. 
The solution exhibits strong recollimation and oscillations. 
It may either reproduce WTTS or CTTS without jet such 
as RY Tau in 2001 when the jet was not visible. Instead, 
a UV shock was dete cted close to the star, possibly cause d 
by the recollimation (jGomez de Castro fc Verdueol l200lt) . 



4. Nonoscillating collimated solution 

4.1. Procedure 

In STT99 and STT02, we have shown that exact MHD 
solutions corresponding to efficient magnetic rotators are 
usually collimated both by the magnetic field and by the 
external gas pressure. However, by decreasing the base pres- 
sure and keeping the other parameters fixed, we obtain a 
limiting solution where the asymptotic structure is colli- 
mated only by the magnetic hoop stress. From underpres- 
sured close to the star, the jet becomes overpressured at 
large distances. For even lower base pressures the wind is 
no longer collimated. We choose in this section a limiting 
solution that gives the widest possible opening of the jet. 
This can be done because the solution is an efficient mag- 
netic rotator with a value of e positive but close to zero 
(cf. STT99). In addition to the parameters defined in Sect. 
13.31 to build our solutions we assumed that the jet speed 
at infinity along the polar axis is Voo = 200 km s~^. 

The process of obtaining a self-consistent solution is a 
long iterating one. As explained in STT02, we start our 
integration at the Alfven point and integrate upstream 
and downstream, each time crossing the critical points and 
changing the parameters until we get as close as possible to 
the guessed values of the parameters. 



4.2. A model for low mass accreting TTauri jets 

Our iterative procedure converges on the following set of 
parameters 

- e = 0.012, 

- S = 0.0778, 

- K = 0.021, 

- iy = 1.5, 

- A = 0.775. 

with a rotation frequency of 

- n^ 5.15 10-6 rad/s, or, V^.o = 8.6 km s'^ 

This last value is in the lower range of the values measured 
for RY Tau. It corresponds to the value of a slowly rotating 
T Tauri star. From this point of view the example of RY 
Tau is an extreme one. 

There is a large discrepancy for RY Tau between 
determinations of the spectroscopic projected rotational 
velocity and of the photometric rotational period. In 
this particular ca se a value o f 52-55 km s~^ is in- 
ferred for Vsini (iPetrov et all . Il999t iMora et all . 1200 it 
I Agra-Amboage et al.l . |2009() . where i is the angle between 
the line of sight and the stellar rotational axis. This value 
of Vsini is among the highest measured for CTTS. This 
specific determination of the projected rotational velocity 
seems to correspond to a special case of obscuration by the 
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Fig. 1. Topolopy of the nonoscillating solution in the 
meridional plane. In a) zoom out showing the collimation. 
In b) zoom close to the star. All distances are given in as- 
tronomical units (AU). The solid (red) lines correspond to 
open fieldlines connected to the star. The dashed (green) 
lines correspond to closed stellar fieldlines and lines con- 
nected to the disk. All open lines are originating in a polar 
cap extending to a colatitude of 15 degrees. The dotted 
circle indicates the Alfven surface. 



disk rather than a real rotational velocity, however. Thus, 
this value may not correspond to the rotation of the star 
itself. Bouvier et al. ()1993lll995[ ) give a photometric period 
of 24 days corresponding to V^^o = 5.1 km s~^, which is a 
lower limit. Flux measurements give a period from 5 to 6 6 
days, closer to typical values for TTS (jPetrov et al.l , 119991 ). 
However, the variation is far from periodic and again one 
cannot exclude veiling from the disk. In other words, ro- 
tation measurements are extremely difficult to obtain for 
this object. We have tried to construct solutions with ro- 
tational values close to the canonical value of 15 kms^^, 
though within a factor of 2 or so, assuming that more pre- 
cise observational data are needed. 

The topology of the streamlines and magnetic fieldlines 
of the solution in the poloidal plane is displayed in Fig. 
[T] The stellar jet is surrounded by a diskwind and is self- 
collimated by its own magnetic field. It is produced in a 
fairly narrow coronal hole of half opening angle 15 degrees. 
The hole is surrounded by a large dead zone extending up to 
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Fig. 2. Plot of the poloidal velocity for the nonoscillat- 
ing solution along the polar axis (solid line) and the last 
streamline connected to the star (dashed) . By construction 
the poloidal velocity is zero on the last streamline where it 
reaches the equatorial plane. This creates the dip seen on 
the logarithmical scale. 



8 stellar radii. This is comparable to the dipolar structure 
of 1.2 kG around BP Tau recon structed from ESPaDOnS 
observations (JDonati et al.l . 120081 ). We also plot in Fig.[2]the 
outflow speed along the polar axis and the one along the 
last open streamline connected to the star. The asymptotic 
speed along the polar axis is 

- Foo = 393 km s^^ , 

which is higher by a factor of 2 than our initial guess. This 
may be quite large, but considering the uncertainties in- 
volved, it can be an acceptable value. 

In Fig. [3] we plot the polar density (a), pressure (b) and 
temperature (c). The pressure is defined up to a constant 
Poi see Eq. (jA.9[) . By varying this constant, we obtain the 
various curves of Figs. ISb and c. Adding different values of 
Po only affects the asymptotic part of the curves. 

Note that the temperature in Fig. [3]: has a maximum 
along the polar axis, reaching a value of one million de- 
grees, which might be a rather high value. However, this 
effective temperature of the plasma is calculated from the 
equation of state of the gas using the total pressure plotted 
in Fig. |3t>. This total pressure may include in addition to 
the kinetic pressu re, ram pressure or Alfven waves, etc. As 
we have shown in lAibeo et al.l ()2007[ ). the effective temper- 
ature can be easily ten times higher than the kinetic one 
with a relative ampli tude of the Alfven waves ^BlB less 
than unity. Note that iGomez de Castro fc Verdrigol (|2007[ ) 
inferred from UV lines high electronic temperatures associ- 
ated with a wind close to lO^K consistent with this scenario 
and an effective temperature of one million degrees. 



For this solution, the physical quantities at the Alfven 
radius are 

- r^ = 9.29ro = 0.104 AU, 

- K = 103knis-i. 

The last open streamline connected to the star corresponds 
to a value of the dimensionless magnetic flux of 

- aout = 0.989. 
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Fig. 3. In (a) we plot of the polar density, in (b) the to- 
tal pressure in dimensionless units along the polar axis, 
including nonkinetic terms such as turbulent magnetic or 
ram pressure. In (c) we plot the polar effective temperature 
for the nonoscillating solution. The three additional asymp- 
totic curves in (c) correspond to three possible values of the 
temperature at infinity, obtained by adding non-zero values 
of the pressure at infinity, as shown in (b). 



We assumed that the observed mass loss rate originates 
completely from the star. It means that the entire mass 
loss rate is contained within a flux tube with its last line 
connected to the star. This assumption yields a mass den- 
sity, particle density and magnetic fleld at the Alfven radius 
of 



g cm 



n^ 


= 3.05 X 10^ cm 


P* 


= 2.48 X 10-15 g 


B* 


= 1.82 G , 
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respectively, where p* is given by (see Eqs. (3.14) and (3.15) 
in ST94) 



M 



wind 



3(5Mwi 



wind 



1 



27rr2T4i/'out ^nrlV, (1 + .Jaout)^/^ - 1 



(12) 



We recall that the mass loss rate of this solution is M^ind — 
3.1 X lO~^M0/yr. From these values at the Alfven radius 
we obtain 

— Poo — 2.72 X lO^^^gcm"'^, deduced from pi, and the 
Alfven Mach number at infinity, 

— rioo ~ 10 cm^"^ assuming the plasma is fully ionized, 

— Bo = 608 G, 

— Boo = 76 mG, 

^ ra7out(f ) = 7.97 To = 0.0891 AU, which corresponds to 
the cylindrical radius where the last line connected to 
the star reaches the equatorial plane (see Fig. [T] b) 

— ^oo,out = 42.8 To = 0.478 AU, 

which corresponds to the radius of the jet at large distances. 
The radii of bow shocks in jet simulations are much larger 
than the jet itself. Th ey can easily be 10 times larger (e.g. 
iMatsakos et al.l . |2009[ ). A jet with a radius of 1 AU may 
create shocks of width 10 AU. Thus the radius we have 
obtained in this solution is consistent with current obser- 
vations, though for low mass accreting stars t he jet radius 
cannot be resolved (|St-Onge fc Bastienl . I2008D . 

The star-braking time calculated from Eq. (jlOp for this 
solution is 



0.6 X lO^yr. 



(13) 



This corresponds to the typical life time of a CTTS 
(jBouvier et all . 119971 ). If k equals 0.2 instead of 0.06, we 
obtain a star braking time of 2 million years. This does not 
drastically affect our final conclusion. This means that the 
jet is efficiently removing angular momentum from the star. 
Thus we may explain the low rotation frequency of those 
particular stars without invoking the disk locking mech- 
anism. This calculation relies on the steady solution we 
used. This does not prevent the star from sudden violent 
outflows caused by reconnection as shown in various numer- 
ical simulations. As we mentioned in the introduction, the 
presence of coronal mass ejections (CMEs) does not affect 
the longterm evolution of the steady ffow. It is interest- 
ing, however , to no te that the braking time obtained in the 
iKiiker et al.l ()2003D simulations is very similar, albeit with 
a different process. The authors conclude that the accretion 
torque spins the star up because the magnetic torque from 
the X-wind is too weak. We propose here alternatively that 
the stellar jet could compensate for it. 

4.3. Higher mass accreting TTauri jets 

In brighter T Tauri jets, higher mass loss rates are usually 
observed. We may extend the same solution in cylindrical 
radius by including the disk- wind component. To obtain 
a jet mass loss rate of Afwind = lO~®M0/yr, the solution 
should be extended to include the fieldline region connected 
to the disk up to adisk.out = 3.07. Note that this value of 
-^wind is in the range of the mass loss rates, deduced fo r 
RY Tau seen in the optical bv lAgra-Amboage et al.l ()2009D . 
typically 0.16 to 2.6 x lO-^Mg/yr. 



The extended jet has a launching region of radius 
T^diak,out — 0.478 AU and an asymptotic cylindrical radius 
of Woo^out — 0.890 AU. The extension of the launching re- 
gion within the disk is quite large. Conversely the expansion 
of the disk outflow is too small. Indeed, we do not expect 
our model to be extended to this distance from the star. It 
does not apply to the Keplerian part of the accretion disk 
region. A more realistic disk wind model would naturally 
give a larger expansion factor of the jet and an asymptotic 
radius of a few tens of AU, as expected observationally. 

We can also keep the same solution considering only 
the stellar wind component but assuming a mass loss rate 
of Mwind = lO~^M0/yr. Accordingly we recalculate the 
physical quantities. Increasing the mass loss rate by a fac- 
tor of 3 just increases the density by a factor of 3 and the 
magnetic field by a factor of \/3- This docs affect the ve- 
locity and the temperature within the jet. In this case the 
values at the Alfven radius are 

- rioo sa 3.210^cm"'^ assuming the plasma is fully ionized, 

- Bo~ 1.1 kG, 

- Boo~ 140 mG. 

This adapts our solution to higher mass loss rates. In these 
cases however, the large observed jet radii imply that there 
is a significant disk wind contribution. Note that this high 
mass loss rate changes the braking time and reduces it by 
a factor of 3 



2 X lO^yr . 



(14) 



It also increases the dipolar magnetic field at the sur- 
face of the star up to 1.1 kG, which is still reasonable 
(jDonati et al.L l2008f) . It does not change the temperature 
profile of Fig. |31 but because the density is significantly 
higher, this would probably lead to observational signatures 
in X-rays. 

5. A recollimating solution 

5.1. Construction of tlie solution 

By changing the initial pressure and the other parameters 
accordingly to keep the same stellar constraints, we obtain 
a second solution, which recollimates at about 38 stellar 
radii. The morphology of this solution is displayed in Fig. 

SI 

The model parameters are close to those corre- 
sponding to RY Tau again and the recollimation dis- 
tance in that case cor responds to the distance of the 
shock seen in the UV by iFerro-Fontan fc Gomez de Castrd 
(I2003D. From the HST o bservations obtained in 1993, 
iGomez de Castro fc Verdu go (2001) deduced the presence 
of a UV shock and suggested that it might correspond to 
the recollimation of the underlying jet. Following this track, 
we suggest that a shock could form in our solution at the 
position where the jet recollimates. Of course, a simple ana- 
lytical model cannot describe the shock itself, this is beyond 
the scope of the present paper. We just note here that the 
outflow velocity at this point is super fast magnetosonic, 
and hence shocks can form without destroying the outflow" . 
This has been shown numerically bv IMatsakos et aD (J2009t) 
for this solution. 

ISt-Onge fc Bastienl ()2008D observed episodic small scale 
outflows in 1998 and 2005 always along the same ejection 
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Fig. 4. Topolopy of the streamlines and magnetic fieldlines 
of the oscillating solution in the meridional plane. In a) 
zoom out showing the recollimation. In b) zoom close to 
the star. All distances are given in astronomical units (AU). 
The solid (red) lines correspond to open fieldlines connected 
to the star. The dashed (green) lines correspond to lines 
connected to the disk. The dead zone here is reduced and 
all lines from the star are opened and not connected to the 
accretion disk. The solution recollimates at 38 stellar radii 
at a distance where the UV shock was observed for RY Tan 
from 1993 to 2001. The dotted circle indicates the Alfven 
surface. 

axis. They proposed that these successive knots merge to 
form the large-scale jet. We suggest then that the two so- 
lutions may correspond to two different stages of the same 
object. The recollimating solution may correspond to the 
stage with a UV shock and the nonoscillating solution to 
another phase where the jet is visible. 



5.2. Physical discussion on the oscillating solution 

After iterating, we found the following set of parameters 

- e = -0.034, 

- (5 = 0.075, 

- K = 0.065, 

- iy = 5.8, 

- A = 0.884. 
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Fig. 5. Plot of the poloidal velocity for the oscillating solu- 
tion along the polar axis (solid line) and the last streamline 
connected to the star (dashed). The jet decelerates after the 
recollimation distance. 



Note that the value of e is now slightly negative. This means 
that the magnetic rotator is important but less efficient in 
collimating. We know from previous study (STT99) that in 
that case the oscillations are more pronounced. 

The corresponding rotation frequency, azimuthal veloc- 
ity and asymptotic speed along the polar axis are 

- f7 = 4.68 X 10-6 rad/s, 

- y^ „ = 7.82 km s-i , 

- Y^^ 186kms-i, 

respectively. This last value is closer this time to our initial 
guess for RY Tau. We plot the polar velocity of this solu- 
tion in Fig. [SI the corresponding polar density, pressure and 
temperature in Fig.[Bl Similarly to Figs.lSb and c, changing 
the value of Po only affects the curves near their minimum 
but not close to the star. 



For this solution, physical quantities at the Alfven ra- 
dius are: 

- r^ = 4.39 7-0 = 0.049 AU 

- V; = 38.75 km s-i . 

We calculate the density again assuming that the observed 
mass loss rate of Mwind = 3.1 x 10~^AfQ/yr entirely comes 
from the star. This means that the jet radius is defined as 
the last flux tube connected to the star that corresponds to 
a magnetic dimensionless flux 

- a{fmtn) = aout = 0.361 (see Fig.H]). 

It corresponds to a mass density, particle density and mag- 
netic field at the Alfven radius of, respectively, 

- rii, — 1.01 X 10^^ cm"'^ is we assume a fully ionized 
plasma as previously , 

- p* = 8.19x 10-14 gcm-3, 

- B^ = 3.94 G . 

From this solution we obtain the following output values 

- Bo = 27.4 G , 

- Boo = 6.28 G, 

- ^ao-out = 2.09 ro = 0.023 AU . 
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Fig. 6. Plot of the polar density (a), the dimensionless 
polar pressure (b) and the polar temperature (c) for the 
oscillating solution. The temperature is the effective one as 
in Fig. [21 The polar effective temperature is fairly constant 
along the jet, suggesting a uniform level of turbulence. The 
three extra asymptotic curves in (c) correspond to three 
possible values of the temperature at infinity, obtained by 
subtracting non-zero values of the pressure at infinity, as 
shown in (b). 



The asymptotic density is higher in this solution because 
we kept the same mass loss rate, but both the velocity and 
the jet radius are smaller. However, the stellar jet radius is 
so small in this case that even if it were drastically enlarged 
by a shock, say at most by a factor of ten, it would remain 
far below any possible detection. 

This solution, as we mentioned, may describe RY Tau or 
any other CTTS during phases where the jet is not visible. 
We suggest that it may also be well adapted to WTTS 
where there is no obvious jet observed despite the similarity 
of the central star. 



Moreover, the star braking time calculated from Eq. 
TOl) is for this solution 



(15) 



T w 7.7x lOV- 



This is longer than the lifetime of CTTS, closer to the time 
spent by the star once the disk is dissipated. Thus the mag- 
netic braking is strongly reduced in this case and inefficient. 
If the accretion rate remains identical, one would expect a 
spin-up of the star. But it is unclear in this solution if the 
accretion rate is not also reduced because of the reduction 
of the large dead zone of the first displayed solution (com- 
pare Fig. [TJj andSJj). 

5.3. A possible application to the RY Tau jet 

By analyzing high spatial resolution images of the clas- 
sical T Tauri star RY Tauri in the T aurus- Auriga cloud 
taken with the Gemini North telescope, ISt-Onge fc BastienI 
(|2008l ) strongly suggested that the RY Tau jet is a pure stel- 
lar jet. They argue that the knots are intermittent microjets 
that merge on longer scales. In 2005, they deduced the tan- 
gential velocity of the brightest knot already seen by HST 
in 1998. Parallely, in 1993 and 2001, a, UV shock w as ob- 
served bv lGomez de Castro fc Verdugol (|2001ll2007t) . They 
interpreted it as a signature of a stellar jet component. 

In view of these observations, it is tempting to suggest 
that the two different MHD solutions discussed above may 
correspond to the two different stages of RY Tau. Thus, to 
explain the succession of knots seen close to the star, the 
first nonoscillating solution would be relevant for this stage. 
During this period, the stellar activity of the central star 
could be enhanced, producing a large visible stellar jet that 
efficiently brakes the central star. 

On the other hand, during periods where the UV lines 
are detected, the outflow could be described by the sec- 
ond oscillating solution with an invisible jet, except for 
the recollimating region where the UV shock is observed. 
Moreover, the width of the stellar jet at the recollimating 
point is around 2 stellar radii, which corresponds to the 
width of the UV zone before the shock, as observed. This 
may take place during a short period, probably linked to a 
lack of magnetic activity, or even during stellar minimum. 

The jet alternately takes the morphology of the two so- 
lutions depending on the activity of the star. A UV shock is 
formed during the low activity phase. Once the activity is 
high, the jet speed increases and blows away the UV shock. 

Note that this multiphase behavior of the large scale 
jet has been found in time-dependent MHD simulations of 
a two-component MH D outflow, when the in itial injection 
speed of the jet varies (jMatsakos et al.l . [20091) . a variabflity 
that may be caused by a magnetic activity cycle. 

The density proflles of ou r two solutions are consistent 
with the densities given by iGomez de Castro &: Verdtigol 
(|2007t ). They measure an electronic density of 10^^ cm~'^ 
in the Silll UV line between 0.006 and 0.3 AU, and 10^° 
cm""^ in the CHI] UV semi-forbidden line between 0.05 and 
1 AU. These two values are consistent with the densities of 
the second solution displayed in Fig.|6^. They also measure 
a density of 10^ cm^'^ in the [OH] forbidden UV line farther 
out between 4 and 100 AU. This last line is not formed in 
the UV shock region. The observed density corresponds to 
the asymptotic density of the flrst solution displayed in Fig. 
111. 



C. Sauty et al.: Nonradial and nonpolytropic astrophysical outflows IX. 



11 



6. Conclusions 

We have used a simple semi-analytical MHD model (ST94) 
to show that the low rotational velocity of T Tauri stars 
can be understood by considering the effect of stellar jets 
to remove their angular momentum. In particular, we mod- 
eled stellar jets from CTTS with a mass loss rate consistent 
with observations. The solutions for coUimated magnetized 
outflows are obtained via a nonlinear separation of the vari- 
ables in the governing full set of the MHD equations. Two 
main results are obtained. 

First, jets can efficiently brake the central star and re- 
move the bulk of its angular momentum over a timescale 
of a million years, at least for low mass accreting stars. 
This corresponds to the typical lifetime of a star in the 
CTTS phase. Furthermore, with this mechanism the star 
can even slow down within 0.6 million years, if the disk 
does not transport angular momentum onto the star. This 
fairly short slow-down time shows the efficiency with which 
a magnetized outflow can remove angular momentum from 
a rapidly rotating young stellar object. This result strongly 
suggests that stellar jets may explain the low rotational 
speeds observed in those objects, independently of any disk 
locking. This mechanism could be complementary to other 
angular momentum removal mechanisms occurring at the 
X- wind emerging from the interaction zone between the 
disk and the magnetosphere. Our solution fairly well re- 
produces the mass loss rate, terminal speed and rotation 
rate of the jet. The effective plasma temperature in the jet 
is reasonable provided that non thermal processes, such as 
Alfven waves, are at work in the plasma of the jets. These 
physical processes are likely to take place in the jets if one 
considers that the magnetic activity of CTTS is comparable 
to, or higher than, that of our own Sun. 

Second, we also analyzed two interesting MHD solu- 
tions. The first corresponds to a gradually cylindrically col- 
limated wide wind without oscillations in its width. The 
second is a narrower outflow that refocuses toward the 
jet axis and also oscillates in width. Note that these two 
very different solutions were obtained by a slight change 
of the parameters. It is a known result in nonlinear equa- 
tions (e.g. Parker solutions, or, their MHD generalization). 
Correspondingly, we discussed the applicability of these two 
solutions in the context of the RY Tau jet, of which re- 
cent observations indicate that the plasma outflow from 
the rotating and magnetized star may have an intermit- 
tency exhibiting at least two phases. During one phase 
the cylindrical wind is relatively wide and nonoscillating. 
During the other phase the outflow is narrow and oscillates. 
Oscillations may produce radiative shocks seen in various 
emission lines, s uch a s the observed UV lines of RY Tau. 
iMatsakos et al.l ()2009() have shown that the self-similar so- 
lutions used in this paper will not be destroyed by the shock 
formation or the disk-wind interaction. 

From a more general perspective, this dichotomic be- 
havior can be seen in the two classes, CTTS and WTTS, 
which may both actually have jets, but these jets are visi- 
ble only CTTS. The connection of CTTS with their disks 
through their magnetospheres may increase the width of 
their jets such as to make them visible. The mass loss 
rate measured in CTTS is thought to decrease as the star 
evolves toward the stage of WTTS. In this later stage, di- 
rect observations are more difficult and upper limits around 
10~^°Mo/yr can only be inferred. Weak T Tauri stars seem 



to have lost the gas component of their accretion disk or at 
least gas falling onto the star. Our second MHD solution 
suggests that even with the same mass loss rate, the jets 
of WTTS may be invisible with the present angular reso- 
lution because of their narrowness. The lack of accretion 
disk prevents the formation of disk winds and the possi- 
bility of higher mass loss rates. The lower jet width may 
be caused either by a low er stellar magnetic activity (see 
IVidotto et al.l . l2009l 120101 ) or by a lack of magnetic con- 
nection into the disk. The magnetic braking time in this 
case would be on the order of 10 million years, which is 
indeed the time a T Tauri star is expected to spend in the 
weak-line stage. 

Our modeling suggests that CTTS may produce visible 
jets in their active phase that would efficiently brake the 
star. Indeed, in the nonoscillating solution the magnetic 
field dominates the dynamics in the jet. Near the star, the 
magnetic field has a dipolar structure and the wind is more 
efficient in this case in extracting stellar angular momen- 
tum. On the other hand, WTTS may have jets as well, 
which are invisible because they are too narrow. These jets 
would be weaker because of the lack of direct connection 
with the accretion disk through a large dead zone. 
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Appendix A: Governing equations for meridional 
self-similar outflows 

The basic equations governing plasma outflows in the 
framework of an ideal MHD treatment for steady, axisym- 
metric flows are the momentum, mass and magnetic flux 
conservation, the frozen-in law for infinite conductivity, and 
the first law of thermodynamics. In particular, the poloidal 
component of the magnetic field can be derived from the 
magnetic flux in spherical coordinates {r,9,(p), 

VA 
B^^—x^. (A.1) 

r smc/ 

We need to specify the latidudinal dependences of the 
velocity, magnetic, density and pressure fields, V, B, P 
and p respectively. This can be resumed in the following 
assumptions (for details see ST94, STT99 and STT02): 
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the square of the poloidal Alfven number, which is a func- 
tion of radial distance only. 

We defined the dimensionless magnetic flux function 
a{R,9) = 2A{r,9)/r'iB^. This quantity is related to G{R) 
through the following expression 
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where G^ {R) is the cross-sectional area of a flux tube per- 
pendicular to the symmetry axis in units of the correspond- 
ing area at the Alfven distance. 

Finally, for homogeneity with the notations in ST94, 
STT99 and STT02, we also introduced the function F(i?), 
which is the logarithmic derivative (with a minus sign) of 
the well know n expansion factor used in solar wind theory 
(|Kopp fc Holzer . 1976) : 



Np, 



l-G^ 



N, 



V 



M 



#2 



G^ D = i-M^ 



(A.19) 



The meaning of the various parameters is discussed in Sec. 
2. 

At the Alfven radius, the slope of AP{R = 1) is p = (2 — 
F^)/t, where r is a solution of the third-degree polynomial: 
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We recall that the value of F defines the shape of the ^""^ ^^^ ^*^^' indicates values at i? = 1 (for details see ST94). 
poloidal streamlines. For F{R) = the streamlines are 
radial, for F{R) > they are deflected to the polar axis 
(F = 2 means cylindrical collimation) , while for F{R) < 
they flare to the equatorial plane. 

Using these assumptions, the usual MHD equation sys- 
tem reduces to the three following ordinary differential 
equations for Il{R), NP{R) and F(R), 
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